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We show that a model where both CP and R-parity are spontaneously broken exists. We study the 
electroweak symmetry breaking sector of the model and find minima consistent with experimentally 
viable Higgs boson masses. We also demonstrate that one can obtain neutrino masses and mixing 
angles within measured values. 
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Breaking CP in the Higgs sector, either spontaneously 
or explicitly, is a theoretically attractive alternative to 
the CP violating phase present in the Standard Model 
(SM) New sources for CP violation from beyond 

the Standard Model scenarios are needed to generate 
enough CP violation to explain the matter-antimatter 
asymmetry in the universe. While in supersymmetric 
models a large number of new phases emerge, in a gen- 
eral minimal supersymmetric standard model (MSSM) 
these phases are strongly constrained by electric dipole 
moments Q. 

A defect of the Standard Model, which persists in the 
MSSM, is that the neutrino masses vanish. Yet the neu- 
trino experiments have provided strong evidence for small 
nonvanishing neutrino mass 0. One popular way to 
explain the neutrino masses is a small violation of R- 
parity 0, iip = ^_i-^3B-l+2s ^ where i3=baryon num- 
ber, L=lepton number, and s=spin of the particle. An- 
other is the seesaw mechanism |^ which can generate 
small masses for neutrinos by allowing Majorana masses 
through the presence of heavy right-handed neutrinos. 

Supersymmetric models share the problem of origin 
of CP violation with the Standard Model. In addition, 
there is no fundamental reason for the existence of R- 
parity in the MSSM, where it is put in by hand in order 
to protect the proton from decaying. However, if only 
lepton number (or baryon number) is violated, the proton 
does not decay. 

Both neutrino masses and CP violation could be ex- 
plained if CP and R were symmetries of the Lagrangian, 
but spontaneously broken by the vacuum. In this letter 
we provide a model where both of these violations are 
intertwined and both are spontaneous. We show that 
while it is non-trivial to satisfy conditions for both sym- 
metry breakdowns at the same time, there are regions in 
the parameter space where we can realize suitable Higgs 
masses, as well as measured neutrino mass differences 
and mixing angles. 

A model where both favored neutrino mass genera- 



tion mechanisms - namely seesaw and R-parity violation 
- with spontaneous R-parity violation was realized in , 
where the spontaneous Rp violation was introduced via a 
term proportional to NLH2. This represents the famil- 
iar bilinear Rp violating term mixing lepton and Higgs 
superfields, LH2, when the right-handed sneutrino field, 
N, develops a VEV. This term also breaks the lepton 
number spontaneously (but not R-parity), and thus in- 
troduces a superfluous massless Goldstone boson into the 
scalar spectrum. The problem can be solved by adding 
a singlet 5" to the theory, through a term N'^S, which 
explicitly breaks lepton number, as S cannot be assigned 
lepton number —2. 

Attempts to violate CP spontaneously, by complex 
VEVs of the neutral scalars exist jlj, but fulfilling the 
experimental constraints has proven difficult. More than 
one Higgs doublet is needed, see e.g. 0. Spontaneous 
breaking of CP is not possible at tree level in the MSSM 
with two Higgs doublets, while it is allowed in a model 
with three doublets Instead of adding doublets, 

one can study extended models, like the NMSSM model 
[lfij |. where the so called /i-problem has been avoided by 
adding a singlet and requiring Z3 symmetry. At tree-level 
one cannot get spontaneous CP violation in this model 
either and consequently radiative corrections were evoked 
[lll |. In that case a very light Higgs boson emerges p^ . 
as also happens in the MSSM, if spontaneous CP vio- 
lation is induced via radiative correction 13]. Another 
possibility studied is to discard the Z3 symmetry com- 
pletely. On one hand, this way one loses the solution to 
the iJ, problem, on the other hand, it is possible to achieve 
SCPV 01 and also solve the problem of domain walls, 
which are created during the EW phase transition as the 
Z3 symmetry is broken spontaneously. 



An interesting model for spontaneous CP violation was 
presented in |l5|. where the Z3 symmetry is replaced 
by R-symmetries on the whole superpotential, including 
non-renormalizable terms [l^ . The method generates 
a tadpole term for the singlet field S in the soft SUSY 
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breaking part of the Lagrangian. This tadpole term al- 
lows for spontaneous CP violation to occur at tree-level 
[lil ll^ . The tadpole is assumed to originate from non- 
renormalizable interactions, which do not spoil quantum 
stability. We adopt this approach here. 
The superpotential of our model is 

W = h]}Q,H2Uj - h%HiQ,Dj - h%HiL,Ej 

where Hi and H2 denote the Higgs doublet superfields, 
Li and Qi the left-handed lepton and quark doublet su- 
perfields, respectively, and Ei and Ui, Di the lepton and 
quark singlet superfields. Right-handed neutrino super- 
fields are denoted by Ni and S is the gauge singlet super- 
field. The terms in the Lagrangian are the only renor- 
malizable ones that respect CP and i?-parity, in addition 
to the gauge symmetry. As all the parameters in the 
Lagrangian are real, this solves the strong CP problem 
The only possible global symmetries are the baryon 
number and Z3. 

The soft SUSY breaking terms in this model are the 
mass terms for scalars and gauginos, and the part mir- 
roring the superpotential with an additional S'-tadpole, 
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+ A'^L,H2Nj + AHH1H2S 
- h.c\ + Ml ''4>14>, + Ma{\l + Af ). 



Here (j) runs over the scalar fields, i,j over the possible 
family indices and a over the three gauge groups. In 
calculations, we take the mass matrices to be flavor diag- 
onal, and assume that only the third generation Yukawa 
couplings (except those for neutrinos) differ from zero. 
We then impose the same texture on the correspond- 
ing ^-terms, and treat the neutrino Yukawa couplings 
as free parameters. The full tree-level scalar potential is 
Vs — ^so/t + Vp + Vd, where Vp and Vd are the usual F 
and D terms. 

The minimization of the scalar potential with respect 
to the fields <j)i and the corresponding phases 6i yield 
constraints later used in finding the scalar mass matrix. 



dVs 



= 0, 



9V, 
39, 



= 0. 



(1) 



Without spontaneous CP violation, the VEVs are real 
and the minimization equations with respect to the 
phases are always satisfied. 

The minimization equations for the charged scalars can 
be trivially solved by setting all charged scalar VEVs 
to zero. As long as the tree-level masses of these fields 
remain positive and the corresponding soft A-terms re- 
main small enough, this is also the global minimum of 



the potential with respect to these fields 0. We use the 
seventeen equations (the phase of the Hi field can al- 
ways be rotated away) given by the moduli and phases 
of the neutral scalar fields to solve for the soft masses 
of the corresponding fields and subset of ^-parameters 
(Aat, As, An-, A^^). Thus the complex VEVs remain 
free parameters and we denote 



(^i> - ( ) ' ^^2^ = ( V2e^''^- ) ' = 



Note that since Rp is violated, the W mass is 
\glv'^ where = vf + vj + « (174 GeV)^. 

Writing the fields as (/) = 0^ + i'Pii with nine neutral 
scalar fields (two Higgs, one singlet and six sneutrinos) 
we get an 18x18 dimensional mass matrix for the scalars. 
The radiative corrections to the scalar masses are imple- 
mented via the one-loop effective scalar potential. 
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where A is the renormalisation scale. A similar term for 
the bottom (s)quark masses contributes significantly, if 
tan (3 is large. Here — ytV2 and to? are the eigenval- 
ues of the stop mass matrix. 
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The loop corrections lead to additional terms in both the 
minimisation conditions and the scalar mass matrix. In 
numerical calculations, we omit the D-term contributions 
and set for simplification, Mq^^ = Mjj^^ = Msusy, with 
MsusY A - 1 TeV. 

The following experimental input is used 

^2 ^ 174 Gev, mw = 80.42 GeV, mf "'"^ = 180 GeV, 
as = 0.102, nir = 1.777 GeV, sin^ 9^ = 0.23124 (4) 

The rest of our free parameters are randomly sampled, 
with sampling ranges as follows (the couplings are con- 
strained by perturbativity): 

0.1 < Ah,7v. < 0.4, 0.2 < As < 0.7, |/i7v| < 10"^ 
0.4 TeV < <^ < 1 TeV, -tt <9^<-k, \{S)\ < I TeV, 
|(i>l)| < lOOkeV, |(iV,)| < 1 TeV, 2 < tan/? < 60, (5) 
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FIG. 1: Lightest physical scalar mass as a function of tan/3. 
The main component of the lightest scalar is Hi 2 (blue rect- 
angle), 5* (light blue diamond), i>L (green triangle), or A'' (red 
triangle) . 



and the A-parameters not eliminated by Eq.JQ) vary be- 



tween < < (1 TeV)/i*^. 



In the limit Aat. — > we recover lepton number conser- 
vation. CP and R-parity are still spontaneously violated, 
and thus in this limit the model contains an experimen- 



There are a total of 55 relevant parameters 
(cr^), e^, As, Xh, Ajv,, h%, M|, As, Ah, An,, A'^ 
and ^. 17 of these can be eliminated by using the min- 
imization conditions, leaving 38 free parameters. Sam- 
pling over such a large space makes it extremely difficult 
to find viable minima of the scalar potential, even with- 
out applying any other cuts or considerations. For the 
full model we choose to set ctr^ ^ = Or^ ^ = to further 
reduce the sampling space. This choice identically solves 
the vacuum conditions dgj^_^ = and = 0. 

In this model, the couplings of the lightest scalar are re- 
duced compared to the SM, due to the (sometimes domi- 
nant) singlet components. Thus the corresponding exper- 
imental limits are smaller as well than the experimental 
lower bound for the SM Higgs boson, mu ^114 GeV, 

see e.g. discussion in [H. In Fig. n we plot the mass of 

the lightest Higgs boson as a function of tan j3, and show 
that experimentally acceptable Higgs bosons masses can 
be found for a large part of the parameter space. 

In a field basis of ,Ni,S, H^, H2, Ao, A3 the neutral 
fermions form the following 11x11 mass matrix: 
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It is easy to see the structure of the usual seesaw mech- 
anism, which produces small neutrino masses m^. 
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mo 



I Mr 



-moMj^^m^, (7) 



where mu <§; M^. Similarly to (7|, there are actually 
several sources for neutrino masses: the usual seesaw and 
the mixing of neutrinos with /i" and gauginos. 

Inspecting the requirement that mo <SC Mr yields 
some qualitative understanding of the model. In partic- 
ular, the left-handed sneutrino VEVs must be small and 
h^iN*) should be of the same order. Thus, although 
(N*) is not bound by any other prior consideration, hav- 
ing /iTv ~ 10~^ results in an upper limit of a few TeV for 
the right-handed sneutrino VEVs. 

We diagonalise M-^o numerically and use Mi ~ M2 ^ 1 
TeV. Great care must be taken, as the elements of M^o 
may vary over ten orders of magnitude, and the eigen- 



values themselves over as much as twenty orders of mag- 
nitude. The diagonalising matrix A/", with AfM^oN'~^ = 
diag(myO , m^. ), has the following general form 



(8) 



Here C) C ^ 1 denote 8x3 matrices that can be deter- 
mined perturbatively, see e.g. 0. Our interest lies in 
the matrix Vi,, the neutrino mixing matrix. Using the 
canonical notation for the neutrino mixing matrix |l9l |. 
we can extract the mixing angles as follows 
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FIG. 2: Correlations of neutrino angles and mass differences. 
Experimentally allowed regions at Scr level are colored. Points 
satisfying neutrino mass constraints are marked as triangles. 
Upper plot: Aniatm vs Attiq, lower plot: sin^ 2023 vs sin^ 613, 
requiring at the same time that the experimental limit for 
sin^ 6112 is satisfied. One point satisfies all constraints (red). 



Once we have sufficient viable minima, i.e. points in 
parameter space that produce a minimum of the poten- 
tial, we apply experimental constraints. We calculate 
how much each minimum deviates from the current ex- 
perimental results and select a number of closest events 
for further analysis. The experimental input we use is 
the following 

sin^ 2023 > 0.89, sin^ 6*13 < 0.047, 



sin^ 612 ~ 0.23 - 0.37, 



712 

atm 
,2 ^ 



"^eV^ - 3.3 X 10" 



Am^ ~ 7.3 X lO^^eV^ - 9.1 x lO^^eV^ 



(10) 



The first stage of finding points for the full model which 
are minima of the scalar potential has a signal-to-noise 
ratio of one over fifty thousand, which shows the odds 
of having an event fall within the constraints. Of these 
points only roughly one in a thousand will satisfy all the 
above constraints. In Fig. |21 we demonstrate that it is 
possible to obtain models where the mass differences and 
experimentally found mixing angles are satisfied. 



Summarizing, we have shown that a model that vio- 
lates both CP invariance and R-parity exists, and we con- 
structed it explicitly. We have shown that experimentally 
viable neutrino and Higgs boson masses can be obtained. 
We have reason to expect that in our model EDM bounds 
and experimental results on kaon (especially ck) and B 
physics can be satisfied j2]| . 
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